submitted to THE ASTROPHYSICAL JOURNAL LETTERS, 07/13/01 



FUNDAMENTAL PHYSICAL PARAMETERS OF COLLIMATED GAMMA-RAY BURST AFTERGLOWS 

A. Panaitescu 

Dept. of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 

AND 

P. Kumar 

Institute for Advanced Study, Olden Lane, Princeton, NJ 08540 



o 
o 

(N 

< 

(N 



> 

in 
o 

oo 

o 

o 

S3 

Or 

6 ■ 

£: 



X 



ABSTRACT 

We determine the basic physical characteristics of eight Gamma-Ray Bursts (GRB) - 980519, 990123, 990510, 
991028, 991216, 000301c, 000926 and 010222 - by modelling the broadband emission of their afterglows. We 
find that the burst kinetic energies after the GRB phase are well clustered around a mean value of 3 x 10 50 
ergs. In contrast, the energy release in 7-rays, after correcting for collimated explosion, varies among bursts by 
more than an order of magnitude. The jet initial apertures are the 2° — 14° range, mildly correlated with the 
energy, half of the jets being narrower than ~ 3°. This implies that, within 100 Mpc, there are about 10 GRB 
remnants (expanding at ~ 0.1 c) which can be resolved with VLBA. For all eight afterglows the total energy in the 
shock-accelerated electrons is close to equipartition with protons. However the slope of the power-law electron 
distribution is not universal, varying between 1.4 and 2.8. In at least half of the cases, the density structure of the 
medium is inconsistent with an r~ 2 profile. A homogeneous medium with density in the 0.1 — 50 cm~ 3 range 
can accommodate the broadband emission of all afterglows, with the exception of 990123, for which we find the 
density to be less than 10~ 2 cm~ 3 . If GRBs arise from the core collapse of massive stars, then such low densities 
indicate the existence of superbubbles created by the supernovae and winds within a cluster of massive stars. 



1. INTRODUCTION 

The detection of very energetic photons (in the GeV range) 
in several Gamma-Ray Bursts (GRBs) and the very short vari- 
ability timescale (Fishman & Meegan 1995), sometimes less 
than 1 ms, exhibited by the 100 keV emission of many bursts 
have lead to the conclusion that they arise from sources that are 
moving at extremely relativistic speeds, with Lorentz factors T 
that could exceed one hundred (Fenimore, Epstein & Ho 1993, 
Piran 1999). Such a highly relativistic motion would follow af- 
ter the release of a large amount of energy in a region of small 
baryonic mass. The non-uniformity in the velocity at which 
various parts of the outflow move leads to internal shocks. In 
this way a fraction of the kinetic energy of the outflow is dis- 
sipated and radiated away in 7-rays (Rees & Meszaros 1994). 
Some of the remaining kinetic energy is converted to electro- 
magnetic radiation when the GRB remnant is decelerated by 
the circum-burst medium (Rees & Meszaros 1992). The result- 
ing external shock energizes the swept-up gas which, similar to 
the GRB phase, emits synchrotron and inverse Compton emis- 
sion, producing an afterglow. Theoretically, the fall-off of the 
afterglow Lorentz factor V is expected to be a power-law in the 
observer time, which leads to a power-law decay of the after- 
glow flux (Meszaros & Rees 1997). This behaviour has indeed 
been observed in many afterglows (e.g. Piran 1999, Piro 2000, 
Wheeler 2000), on timescales of days. 

Due to the relativistic beaming of the emission, the observer 
receives radiation mostly from the fluid moving within an angle 
1 /r radians off the observer's line of sight toward the fireball 
center. Thus as T decreases, the size of the "visible" region 
increases. If the GRB remnant is collimated into a jet, then at 
some time tj the entire jet surface becomes visible to the ob- 
server. This time is given by 



tj = 0.4 (z + 1) (£0,50^) V3 0o,-i d ay 



(1) 



z being the burst redshift, £?o,50 the initial jet energy measured 
in 10 50 erg, #o.-i its initial half-opening measured in 0.1 radi- 
ans, and no the external medium density in cm~ 3 . The lack of 
emitting fluid outside the jet opening leads to a faster afterglow 
decay after tj. Furthermore, around tj the jet starts to expand 
laterally (Rhoads 1999), its sweeping area increases faster than 
before, leading to a stronger deceleration and increasing even 
more the afterglow dimming rate. The most important signa- 
ture produced at tj by the jet collimation is the achromaticity 
of the afterglow emission break, which distinguishes it from the 
chromatic light-curve steepening that the passage of a spectral 
break through the observing band would yield. So far, there are 
eight well observed GRB afterglows (980519, 990123, 990510, 
991208, 991216, 000301c, 000926, and 010222) for which a 
break or a steep decay has been observed in their optical light- 
curves. These are the afterglows we model in this paper. 

The afterglow flux at a given frequency and time depends on 
the jet speed, on the properties of the external medium, and on 
the micro-physics of relativistic shocks. Our aim is to constrain 
some of these properties by modelling the afterglow dynamics, 
calculating its emission and using the observational data to de- 
termine i) the jet energy and collimation, ii) the efficiency at 
which relativistic shock accelerate electrons and generate mag- 
netic fields, and Hi) the external medium type and density. 

2. THE AFTERGLOW MODEL 

There are some simplifying assumptions made in our numer- 
ical modelling, which allow the high computational speed nec- 
essary for a parameter space search. The most important sim- 
plification is that the jet front is homogeneous: the energy per 
solid angle and Lorentz factor have the same value in any direc- 
tion within the jet aperture, and the internal energy density of 
the shocked fluid is uniform, at the value set by the shock jump 
conditions. Implicitly, we also assume that the jet has a sharp 
edge. 
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The dynamics of the jet is calculated by tracking its energy 
(some of which is lost radiatively), mass and aperture, which 
increase as the jet expands and sweeps-up the surrounding gas 
(Kumar & Panaitescu 2000, Panaitescu & Kumar 2000). There 
are only three parameters that give the jet dynamics - the initial 
jet energy Eq, initial half-angle 9q, and external particle density, 
n for a homogeneous medium or the constant A for a wind-like 
medium with profile n(r) = Ar~ 2 . 

The calculation of the afterglow synchrotron and inverse 
Compton emission is based on the assumption that the electron 
distribution M e (j) produced by the shock acceleration in the 
downstream region is a power-law N e oc 7~ p (this is supported 
by the observed power-law decay of many afterglows), starting 
from a minimum electron energy 7im e c 2 and ending at a high 
energy break 7*m e c 2 . The existence of this high energy break 
is due to the escape of particles from the acceleration region, to 
radiative losses, and is also required to have a finite energy in 
electrons if p < 2. For simplicity we approximate the cut-off at 
7* as a steeper power-law of index q. 

There are three or five model parameters pertaining to the mi- 
crophysics of shock acceleration and magnetic field generation, 
based on which we calculate the co-moving frame afterglow 
spectrum. The magnetic field strength is simply parameterized 
by the fraction eb of the post-shock energy density that resides 
in it. The distribution of the injected electrons is defined by the 
fractional energy e e in electrons if they all had the same 7 = 73 
(thus e e parameterizes 7,) and the index p of the power-law 
distribution above 74. The cut-off 7, is set by the the fractional 
energy e of the electrons between 7; and 7, . 

The synchrotron spectrum is piece-wise power-law, with 
breaks at the self-absorption frequency v a , injection frequency 
v. L corresponding to the minimum electron 7,-, cooling fre- 
quency v c corresponding to an electron cooling timescale equal 
to the dynamical time, and cut-off frequency z/* associated with 
7* . The afterglow spectrum and light-curve is determined by 
the evolution of these spectral breaks and of the flux at peak 
frequency v v — min(V c , Vi). We note that for a spreading jet 
(t > tj e t), the afterglow emission at frequencies above z/j de- 
cays approximately as F v oc t~ p , irrespective of the location 
of v c , provided that the electron cooling is dominated by syn- 
chrotron emission. Therefore the temporal slope of the post jet- 
break optical or radio light-curves gives a direct measurement 
of the electron distribution index p. 

We calculate numerically the afterglow emission for an ob- 
server lying on the jet axis, by integrating the synchrotron and 
inverse Compton emissions over the jet dynamics, taking into 
account the variation of the bulk Lorentz factor T on the equal 
arrival time surface. We note that the afterglow light-curve has a 
weak dependence on observer location offsets less than 60 ■ The 
equations we employ for the calculation of the jet dynamics and 
of the received radiation (Kumar & Panaitescu 2000, Panaitescu 
& Kumar 2000, 2001) are valid in any relativistic regime. For 
the parameters we find for various GRB afterglows, the non- 
relativistic regime should set in after the last available obser- 
vations. It is the mildly relativistic regime that corresponds to 
observations made later than about 10 days after the GRB event, 
for which a numerical treatment is most appropriate. 



3. FEATURES OF GRB JETS 

Based on the afterglow model outlined above, we model the 
broadband, time-dependent emission of eight GRB afterglows: 
980519, 990123, 990510, 991208, 991216, 000301c, 000926 
and 010222, whose light-curves exhibited a break, offering thus 
the possibility of determining the initial jet opening 8q and en- 
ergy -Eo-Q The basic afterglow parameters: Eq, 6q, n, e e , Eb, 
and p (as well as e and q when they are relevant) are determined 
by minimization of \ 2 (maximization of the likelihood to ob- 
tain the observed fluxes). Figure 1 shows the data and the model 
light-curves in three distinct domains: radio, optical and X-ray. 
The averages and dispersions of parameters for our eight GRBs, 
calculated from the best fit parameters of individual bursts, are 



E = (2.6 ± 1.1) x 10 50 erg 
log(n/cm- 3 ) = 0.13 ± 1.35 , 
log£ fl = -2.4±1.2, 



P 



= 6.1° ±4.5° , 
= 0.062 ±0.045 
1.87 ±0.51 . 



Figure 2 shows the best fit parameters and their 90% con- 
fidence levels, determined by variation of \ 2 around its min- 
imum, for our set of eight afterglows. As illustrated, we find 
that jet energies that are clustered between 10 50 and ~6x 10 50 
erg,[] somewhat smaller than the typical kinetic energy of a su- 
pernova. Thus only a very small fraction (~ 10~ 4 ) of the en- 
ergy budget expected for GRB progenitors (Paczyhski 1998, 
Meszaros, Rees & Wijers 1999, MacFadyen, Woosley & Heger 
2001) is given to highly relativistic ejecta. It is then surprising 
that the variation of this fraction among bursts is merely a factor 
of a few. 

We have found that models with a homogeneous medium can 
explain well the broadband emission of all eight afterglows, al- 
though in one case (000926) the best fit is unsatisfactory (% 2 
is rather large). For 980519 and 990510 we determine external 
particle densities around ~ 0.1 cm~ 3 , typical for warm inter- 
stellar medium, while for 991208, 991216, 000301c, 000926 
and 010222 the densities we find are around 10 cm -3 . For 
990123 we obtain an external density below 10~ 2 cm~ 3 , simi- 
lar to the afterglow 980703 (Panaitescu & Kumar 2001). 

Our afterglow calculations also show that an r -2 density pro- 
file of the external medium is compatible with the emission of 
991208 and 991216, and can marginally accommodate the af- 
terglows 000301c and 010222, but cannot explain the broad- 
band emission of 980519, 990123, 990510 and 000926. There- 
fore our findings are not consistent with the "pure" winds ex- 
pected from Wolf-Rayet stars (Chevalier & Li 1999) and are in 
accord with the results of Ramirez-Ruiz et al. (2001), who have 
shown that the interaction between such winds and the circum- 
stellar medium leads to the formation of a quasi-homogeneous 



shell extending from ~ 10 cm to > 10 cm 

Regarding the microphysics of shocks, the results presented 
in Figure 2 indicate that i) the total electron energy is close 
to equipartition with protons, it) the magnetic field strength, 

1/2 

B oc Eg , is within two orders of magnitude of the equipar- 
tition value, and Hi) the power-law distribution of shock- 
accelerated electrons does not have an universal index p, with 
1.4 < p < 2.8. Given that the post jet-break emission falls-off 
as ir p , the shallow decays of the optical emission of 010222 
and radio light-curves of 991208, 991216 and 000301c require 

1 Other afterglows did not have such a break for the duration spanned by the observations. Although they may have been collimated, the initial jet angle cannot be 
determined without having an observational constraint on the jet break time tj . 

2 Using the analytical expression for the jet break time, Frail et al. (2000) have calculated jet apertures 9o for more afterglows and found that the energies released 
in 7-rays by GRB jets are also clustered within a factor of 4. 
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models with p ~ 1.5. Consistency between such hard electron 
distributions and the observed optical spectral slopes implies 
that, for these four afterglows, the cooling frequency v c was 
below the optical domain. 

In our modelling, the passage of a high energy spectral break 
v* yields most of the steepening of the optical decay observed 
in the afterglows 991208, 991216 and 000301c. This allows 
us to determine the fractional electron energy up to 7* corre- 
sponding to v*. We find this fraction to be in the range 1/3-2/3 
(Figure 2), close to equipartition, which provides a natural rea- 
son for the existence of the break at 7* . We note that hard elec- 
tron distributions and equipartition electron energies find mu- 
tual consistency in the shock acceleration treatment of Malkov 
(1999). 

Assuming an uniform jet, the jet energies and apertures in- 
ferred here and the observed 7-ray fluences require an effi- 
ciency of the 7-ray mechanism in dissipating the jet kinetic 
energy and radiating it in the 20 keV-1 MeV range which ex- 
ceeds ~ 50%. However if the 7-ray emission arises from bright 
patches caused by angular fluctuations (Kumar & Piran 2000) 
of the kinetic energy on the jet surface, the actual efficiency 
can be much smaller, closer to that obtained from numerical 
modelling of internal shocks (Spada, Panaitescu & Meszaros 
2000). These bright spots could also induce short timescale 
fluctuations in the early X-ray and optical afterglow emission, 
before these fluctuations disperse, and may have already been 
observed in the optical light-curve of 000301c. 

4. CONCLUSIONS 

The findings reported above should shed some light on the 
GRB progenitors. Homogeneous media of low density, below 
10" 2 cm- 3 , as found for the afterglows 990123 and 980703, 
indicate a galactic halo or a hot component of the interstel- 
lar medium. External particle densities in the 0.1 — 1 cm" 3 
range, as found for 980519 and 990510, are characteristic for 
the interstellar medium, while larger values of ~ 10 cm -3 , as 
determined for the afterglows 991208, 000301c, 000926, and 
010222, are typical for diffuse hydrogen clouds. Even larger 
particle densities, above 100 cm" 3 , would be consistent with 
the undisturbed, dense molecular clouds expected in the cur- 
rently popular collapsar model (Woosley 1993, Paczynski 1998, 
MacFadyen & Woosley 1999) for GRB progenitors. Scalo & 
Wheeler (2001) have emphasized that the winds and super- 
novae occurring in a cluster of massive stars create 10 pc-1 
kpc superbubbles whose density can be as low as 10" 3 cm" 3 . 



Furthermore they argued that variations in the cluster age and 
the density of the giant molecular cloud into which the super- 
bubble expands, as well as the interaction with the winds from 
other clusters, may yield circumburst medium densities span- 
ning a few orders of magnitude. 

The relativistic kinetic energies of the eight GRBs analyzed 
here show a remarkably narrow distribution, with mean value of 
~ 3 x 10 50 erg and dispersion of ~ 10 50 erg. On the other hand, 
the distributions of the 7-ray energy output (assuming uniform 
jets) and of the jet initial aperture are significantly broader, with 
a width of one order of magnitude. 

That half of the afterglows analyzed here are narrower than 
~ 3° (Figure 2), indicates that the distribution of initial jet an- 
gles is dominated by very narrow jets.^ Taking into account that 
star formation rate at redshift z ~ 1 was an order of magnitude 
larger than at present (Madau, Pozzetti & Dickinson 1998), the 
rate of collimated GRBs is <; 10 3 times smaller than that of su- 
pernovae. As pointed out by Paczynski (2001), if the fraction 
of core collapse supernovae associated with GRBs is at z ~ 1 
the same as in the nearby Universe then, within 100 Mpc, there 
should be < 10 GRB remnants that are sufficiently bright and 
large after several years to allow the resolution of their non- 
spherical structure with VLBA. These remnants will have typ- 
ical velocities of 0.1 c, therefore the expansion of the nearest 
remnants, over several years, could be detectable. 

From the afterglow parameters presented in Figure 2, one can 
show that the jet Lorentz factor at the end of the GRB phase is 
160 ± 70. For the typical jet energy given above, this implies 
that the baryonic mass ejected at ultra-relativistic speeds by the 
GRB explosion is about 1O" 6 M0. If the GRB progenitor is a 
10 Mq star, the largest jet apertures (~ 10°) and the baryonic 
loads determined here imply that the region of the star through 
which the jet propagates is evacuated up to 1 part in 10 5 . 

Our calculations led the surprising result that the electron en- 
ergy distribution index varies significantly among GRB after- 
glows, suggesting that the details of the acceleration of particles 
at ultra-relativistic shocks are not universal. We also find that 
the energy imparted to electrons at relativistic shocks should be 
close to equipartition with protons. 
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FIG. 1.— Radio, optical and X-ray emission and model light-curves for the GRB afterglows 980519, 990123, 990510, 991208, 991216, 000301c, 000926, and 
010222 (legend of middle graph applies to all panels). The numerical light-curves have been obtained by minimization of \ 2 between model emission and the radio, 
millimeter, sub-millimeter, near infrared, optical, and X-ray data (only a part of the used data is shown in this figure). The parameters of each model are given in 
Figure 2. Optical data has been corrected for Galactic dust extinction. The spread around the model curves exhibited by the radio emission of 980519, 991208, 
991216, 000301c and 000926 can be explained by fluctuations due to scatterings by the inhomogeneities in the Galactic interstellar medium (Goodman 1997). Fluxes 
have been multiplied by the indicated factors, for clarity. 
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FIG. 2. — Best fit model parameters and their 90% confidence intervals for the afterglows shown in Figure I. a: jet initial energy Eo, b: initial half-angle 8o, 
c: external medium density n, d: parameter e e for minimum injected electron energy (filled circles), fractional electron energy e up to the 7* -break (filled squares, 



triangles indicating lower limits), and the total electron energy - — = e e if the electron distribution extends to infinity (open squares), e: magnetic field parameter eg, 

P z 

f: index p of the power-law injected electron distribution (filled circles), and index q above the 7, -break (filled squares). For 010222 we find an unusually small 
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